A Carbon Nanotube Based Nanorelay 



(N 
O 

o 

(N 

wo: 
< 

(N 
(N 



i 

-a 
c 

o 
o 



> 

00 

o 

O 

-i— > 
c3 



i 

C 

O 

o 



X 



J. M. Kinaret, T. Nord and S. Viefers 
Department of Applied Physics, Chalmers University of Technology 
and Goteborg University, SE-^12 96 Goteborg, Sweden 
(February 1, 2008) 

We investigate the operational characteristics of a nanorelay based on a conducting carbon nan- 
otube placed on a terrace in a silicon substrate. The nanorelay is a three terminal device that acts 
as a switch in the GHz regime. Potential applications include logic devices, memory elements, pulse 
generators, and current or voltage amplifiers. 



Nanoelectromechanical systems (NEMS) are a rapidly 
growing research field with substantial potential for fu- 
ture applications. The basic operating principle under- 
lying NEMS is the strong electromechanical coupling in 
nanometer-size electronic devices in which the Coulomb 
forces associated with device operation are comparable 
with the chemical forces that hold the devices together. 
Carbon nanotubes (CNT)EJ are ideal candidates for nano- 
electromechanical devices due to their well-characterized 
chemical and physical structures, low masses, exceptional 
directional stiffness, and good reproducibility. Nanotube- 
based NEMS have internal operating frequencies in the 
gigahertz range, which makes them attractive for a num- 
ber of applications. Recent progress in this- direction in- 
cludes fabrication of CNT nanotweezers,0a CNT bas< 
random access memoryjj and super-sensitive sensors.B' 

In this paper we consider another example of CNT 
based NEMS, a so-called nanorelay. This three-terminal 
device consists of a conducting CNT placed on a terraced 
Si substrate and connected to a fixed source electrode. A 
gate electrode is positioned underneath the CNT so that 
charge can be induced in the CNT by applying a gate 
voltage. The resulting capacitive force between the CNT 
and the gate bends the tube and brings the tube end 
into contact with a drain electrode on the lower terrace, 
thereby closing an electric circuit. We describe the, sys- 
tem with a model based on classical elasticitgJjheoryu and 
the orthodox theory of Coulomb blockade,tm and study 
its IV-characteristics and switching dynamics. Theoret- 
ical studies of a related— two-terminal structure have re- 
cently been reported.E90 

Model system. The geometry of the nanorelay is de- 
picted in Fig. |l|. We model the CNT as-, an elastic 
cantilever using continuum elasticity theory:0 Assuming 
that only the lowest vibrational cigenmode is excited, 
and that the bending profile upon applying an external 
force is the same as that of free oscillations, one can 
express the potential energy of the bent tube in terms 
of the deflection x of its tip as V = kx 2 /2. The ef- 
fective spring constant k depends on the geometry of 
the tube and is approximately given by k rs 3EI/L 3 . 
Here E is Young's modulus, (paaepmcntally determined 
to be approximately 1.2TPa,t3'ta L is the tube length 
and / = tt(D* — D|)/64 its moment of inertia, D Q and 
Di being the outer and inner diameters of the (multi- 



walled) CNT. The effective mass of the tube is m e // = 
k/fl 2 w 3M/(1.875) 2 , where M. is the total tube mass 
and its lowest eigenfrequency.13 It is known experimen- 
tally that-Q-factors of CNT cantilevers are of the order of 
170-500£fl We model this by a phenomenological damp- 
ing force — jdi in the equations of motion. 




FIG. 1. Schematic picture of the model system consist- 
ing of a conducting carbon nanotube placed on a terraced 
Si substrate. The terrace height is labeled h, and q denotes 
the excess charge on the tube. The CNT is connected to 
a source electrode (S), and the gate (G) and drain (D) elec- 
trodes are placed on the substrate beneath the CNT at lengths 
L and Lc away from the terrace. The displacement x of the 
nanotube tip is measured towards the substrate. Typically, 
L « 50 - 100 nm, h ps 5 nm. 



For a metallic nanotube, the effective impedance Z is 
dominated by the contact resistance at the source con- 
tact and is mostly ohmic. If a gate voltage is applied, 
the resulting capacitive force bends the tube so that a 
tunneling current can flow to drain, whereas due to the 
long tube-gate distance no electrons tunnel between the 
tube and the gate electrode. Finite element calculations 
show that the drain- and gate capacitances are well ap- 
proximated by those of parallel plate capacitorsJla We 
thus model them as 



r m - 2Co 

WW - i_ K *(i_c /C h )- 



(i) 



Here Co = Cd(x — 0) is the drain capacitance for a 
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horizontal tube and Ct = Cg(x — h) is that for the 
tube in contact with the drain eleotaode. The latter 
can be estimated from experiments. 113 The constant k 
(0 < k < 1) accounts for the fact that the deflection 
at the position of the gate is smaller than that of the 
tip. The corresponding capacitive forces take the form 
F c = -(Q 2 /2)V(1/C(xj). 

The tunneling resistance is modeled as Rt = 
Ro exp[(/i — a;)/A], where the tunneling length A depends 
on the contact material and is typically of the order pf 
0.5 A, and Rq is estimated from experimental results.E3 
Following Ingold and NazarovO for the case of zero tem- 
perature and an ohmic environmental impedance Z, we 
determine the tunneling rate at a given deflection x and 
source-drain voltage V s d as 



T(x) 
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where P{E) 1 the probability for energy exchange between 
the tunneling electron and the environment, is deter- 
mined self-consistently from 



EP(E) = - / dE' 
9 Jo 
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Here, E c (x) = e 2 /(2C(xj) and g = R K /Z where R K = 
h/e 2 25.8 kO is the von Klitzing constant. 

The above model results in a set of coupled non-linear 
differential equations describing both the mechanical mo- 
tion and the current flow in the system, 



m eff x(t) 
Zq(t) 
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where it is the stochastic tunneling current, and the zero 
of the potential has been chosen such that Vd — 0. The 
contact force F contac t represents the short range forces 
between the tube tip and the drain electrode important 
at large deflections; presently, we model the surface in- 
teractions as elastic collisions. For a set of applied volt- 
ages, design parameters and initial conditions, we solve 
the equations of motion numerically using a fourth or- 
der Runge-Kutta method. The tunneling probability is 
computed at every time step, and tunneling events are 
treated as a (weighted) random process (for algorithmic 
details, see Ref. . 

Results. In our simulations we have studied the dc IV- 
characteristics of the nanorelay and its response I d (t) to 
a gate voltage pulse. Fig. || shows the drain current as a 
function of gate- and source-drain voltages. Due to the 
exponential dependence of the tunneling resistance on 
tube deflection, there is a sharp transition from a non- 
conducting (OFF) to a conducting (ON) state when the 
gate voltage is varied at fixed source-drain voltage. The 



sharp switching curve allows for amplification of weak sig- 
nals superimposed on the gate voltage; detailed analysis 
of the amplification characteristics will be pursued else- 
where. The inset shows a high resolution plot of the drain 
current as a function of V sg for various V sd , and demon- 
strates that the pull-in voltage V* depends weakly on 
V sd as expected, giving rise to highly non-linear Id(V s d) 
characteristics near the transition. 




FIG. 2. Drain current as a function of source-drain voltage 
and gate voltage. The parameters used here are L = 75 nm, 
h = 3nm, Di — 2nm, D = 3nm, Lq = 0.7 L and Z — 8kf2. 
The pull-in voltage V* g for various values of V a d can be read 
off the higher resolution IV curves shown in the inset. 



For certain choices of parameters, hysteresis occurs in 
the OFF-ON-OFF transition. This is illustrated in Fig. | 
which shows I d as a function of V sg for up- and down- 
sweeps. The hysteresis is due to the appearance of two 
stable cantilever positions for certain design parameters 
and gate voltages: one of the stable positions is at con- 
tact (ON state) and the other is a slightly bent tube 
configuration (x < h/2) corresponding to an OFF state. 
The inset of Fig. || shows the contours of zero net force 
as a function of gate voltage and tube deflection (net 
force is positive to the right of the contour) for two sets 
of design parameters; hysteresis arises for the first set of 
parameters in a range of gate voltages when there are 
two zero force configurations. The amount of hysteresis 
can be controlled by device design — for example, plac- 
ing the gate electrode further from the drain decreases 
hysteresis. The hysteretic switching characteristics may 
be utilized to construct nanotube-based nanoelectrome- 
chanical memory elements. 

Finally, we studied the switching dynamics of the 
nanorelay, choosing a set of parameters for which hys- 
teresis is unimportant. The simulation presented here 
started in the OFF state with V sg = 0V, switched to the 
ON state (V sg — 3V, see Fig. ||) at t = 0.5 fis and back to 
V sg = 0V at t = 1 /is. The resulting switching times were 
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potential) and vdW forces (based on a Lennard-Jones po- 
tential) to the mechanical and electrostatic ones, and find 
that the sticking problem can be avoided by a suitable 
choice of parameters such as shorter or thicker CNTs. 
Another effect we have neglected is the possibility of field 
emission,! 2 !! which may play a role in the geometry con- 
sidered in this paper. An investigation of these points 
will be the focus of a future publication. 

In conclusion, we have performed a theoretical analy- 
sis of a CNT-based three-terminal nanoelectromechanical 
switch. The large parameter space for its design allows 
for a number of potential applications, such as amplifiers, 
logic devices or memory elements. 



FIG. 3. Switching curves for up- and down-sweeps of V sg 
for various V s d showing hysteresis (L = 75 nm, Lg = 0.7L, 
Di — 2nm, D — 3nm). Inset: contours of zero net force as 
a function of V 3g and cantilever deflection x for two devices 
(total force is positive to the right of each curve). The ap- 
pearance of two configurations with zero net force for some 
Vsg gives rise to two stable cantilever positions corresponding 
to ON and OFF states, with associated hysteresis. Device 
showing hysteresis L = 75 nm, Lg = 0.7L, device without 
hysteresis L — 100 nm, Lg = 0.5L, other parameters as given 
above. 
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tens of nanoseconds for the OFF-ON transition (e.g. rise 
time tr w 20ns), while the ON-OFF transition was con- 
siderably faster (fall time rp « 0.1ns). This asymmetry 
is in part due to the large difference between the step 
height and the tunneling length, in part to the cantilever 
bouncing off the drain surface during the OFF-ON tran- 
sition. The switching dynamics are quite sensitive to the 
geometrical parameters which means that they can be 
optimized for specific applications. The time constants 
also depend on the details of the surface force F contact 
which are beyond our model. In particular, the OFF-ON 
transition time is greatly reduced by dissipative surface 
processes [e.g. phonon emission) that dampen tube vi- 
brations, and allow the CNT to settle faster at a station- 
ary near-contact position, ._. ._. 

The van der Waals (vdW)Efl and adhesiveO forces be- 
tween the CNT and the_saibstrate were neglected in our 
model. Dequesnes et a/£d recently studied the effect of 
vdW forces on the pull-in voltage V* for a two-terminal 
device (the entire substrate acting as a gate electrode). 
The main result is that vdW forces reduce the pull-in 
voltage, but do not change the qualitative features of the 
OFF-ON transition. Moreover, the significance of the 
vdW forces decreases with decreasing tube length and 
increasing terrace height or tube diameter, and their ef- 
fect on V* is negligible for a wide range of design pa- 
rameters. Concerning the ON-OFF transition, adhesive 
forcestJ, enhance the tendency of the CNT to remain in 
contact even after the gate voltage has been switched 
off. We have performed a simple analysis comparing the 
magnitude of the adhesive forces (modeled by a Morse 
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